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Abstract 

O 

^^ I Numerical simulations of phase ordering under dissipative dynamics in a (2+1)- 

^SJ I dimensional 3-vector model with 0(3) symmetry are reported. The energy func- 

pi i' tional includes terms which stabilize the size of extended topological defects. 

They emerge at the end of the coarsening process as particle- or antiparticle- 
like structures floating in the globally aligned vacuum. Approximate power-law 
growth of disoriented domains (with an exponent near 0.4) is found to be rather 
insensitive to the size of the defects. An optional filter for conservation of wind- 
ing number allows to study phase ordering in high defect-density environment 
^O ■ which leads to large clusters of particles surrounded by the aligned field. 
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1. Introduction 

Effective field theories for order-parameter A^- vector fields can be powerful tools for the study 
of the dynamical behaviour of complex microscopic many-body systems. Standard examples 
are the local magnetization 3- vector field $ defined over (2-1-1 )-dimensional space-time for the 
description of pseudo-2-dimensional spin systems in solid state physics [|ll; similarly, 4- vector 
fields defined over (3-|-l)-dimensional space-time model low-energy QCD in chiral effective field 
theories for the dynamics of pions and their interactions with gauge fields [|2| . Extending these 
models to finite temperatures T may allow to investigate features of ordering transitions in 
spin systems or, for the hadronic case, in the cooling phase of the early universe or immediately 
after a heavy-ion collision. 

It is expected that during such a cooling process the system undergoes a transition from 
a hot state in which the global symmetry of the effective action is manifest into a cold state 
with spontaneously broken symmetry, i.e. with global alignment of the ordering vector field 
in a randomly chosen direction. The randomly oriented aligned domains which characterize 
transient intermediate stages of this ordering transition have recently found increased attention 
for the case of the chiral meson field, because it was suggested that they might cause anomalies 
in the multiplicities of emitted pions which could serve as a signature for the phase transition 
itself 0. Similarly, the density of topological defects in the field configurations has been related 
to production rates of extended particle and antiparticle structures embedded in the aligning 
field i. 

This latter aspect is based on the particularly interesting topological properties of (2-|-l)D- 
0(3) and (3-|-l)D-0(4) models. For appropriate boundary conditions the field configurations 
fall into separate classes which can be characterized by integer topological winding number B. 
For T = the classical ground state configuration (the 'vacuum') and the structure of the 
particle-like excitations (the 'solitons' or 'baryons') are obtained through minimization of the 
corresponding static energy functional -E[$]. The vacuum consists of a globally aligned field 
with constant length $ = /q and B = 0. In case of exact 0{N) symmetry different global 
orientations are degenerate. Explicitly symmetry-breaking terms or boundary conditions select 
a specific orientation as the true vacuum. Static topological defects are local minima of E[^] 
in sectors with B ^ 0. They represent ordered configurations with large gradients in direction 
$ and length $ of the field vector. Large deviations of $ from the vacuum value /o are denoted 
as 'bags'. Stabilizing terms present in E[^] determine soliton and bag size, shape and spatial 
profile of isolated defects. 

Such localized extended configurations have been identified with particle-like structures, 
low- lying excitations with smoothly varying magnetization in Quantum Hall Ferromagnets |^, 
or baryons embedded in and interacting with the chiral mesonic field Q. In these cases the 
corresponding local winding density has been interpreted as local electric charge density, or 
baryon density, respectively. 

Topologically nontrivial configurations can unwind at space-time points where the length 
$ of the field vector vanishes. So configurations (with suitable boundary conditions) can be 
specified by a definite winding number B only if there is no point at which $ = 0. Representing 
local energy minima, static B ^ configurations are separated from global vacuum by barriers; 
they may be destabilized by symmetry breakers which reduce these barriers. In time-dependent 
(nonequilibrium) processes even for exact 0{N) symmetry B might change through fiuctuations 



or evolutions where $ passes through zero at some point. To study the dynamical behaviour 
of field configurations at nonzero temperature for given average soliton density and fiuctuating 
B therefore would require the inclusion of a chemical potential and the construction of a large 
grandcanonical ensemble with numerous individual events in order to allow for a well-defined 
average 5-number. Alternatively, however, if different B sectors are separated by barriers, 
we can put a constraint on B. This is very naturally implemented in lattice simulations if in 
a time evolution the updates with occasional jumps in B are rejected. Such a 5-filter then 
allows to investigate phase ordering in high defect-density environment without the need for a 
grandcanonical ensemble. 

The time evolution of classical fields on a lattice is considered as the evolution of low- 
frequency modes subject to the noise of ehminated high frequency fiuctuations. In addition 
to stochastic forces this provides a dissipative (first-order) time-derivative term which removes 
energy from the field configurations and drives the system to the thermal equilibrium at the 
temperature set by the noise term. This temperature itself may change in time as a consequence 
of the spatial evolution of the system, or as externally imposed quench. Numerical simulations 
then may follow individual 'events' with fixed winding number B which start out from an initial 
random configuration (which would represent a member of an ensemble of high temperature) 
with correlation length less than the lattice constant. In such initial configurations the field 
vectors vary randomly from one lattice point to the next, therefore these configurations are 
characterized by large local gradients, i.e. large (positive and negative) local winding density, 
and large total energy. 

Ordering proceeds through growth of domains which comprise in their interior increasing 
numbers of lattice vertices with aligned field vectors. In the absence of any explicit symmetry 
breaking the relative orientation of the aligned field for different domains is random. Therefore 
large field gradients then are confined to the boundaries of these 'disoriented' domains. The 
soliton stabilizing terms in £'[$] then lead to a dynamical interplay between the growth of 
disoriented domains and the formation of energetically favorable localized extended structures. 
Finally, we expect a few well-developed solitons to remain embedded in an otherwise fully 
aligned vacuum, each soliton with its own (approximately integer) winding number Bi which 
add up to the conserved total integer winding number B. Residual interactions and (at finite 
temperatures) remaining field fiuctuations cause slow motion of the localized structures such 
that on very long time scales some of them may meet and combine or annihilate to form larger 
or smaller B.^ + Bj structures. 

Evidently, the whole evolution of one individual configuration is an extremely complex 
process, and after all, only statistical statements averaged over many events will be of interest. 
However, for a sufficiently large lattice which finally still contains a large number of individual 
solitons, lattice averages for individual events will already provide good approximations to 
ensemble averages. In the present work we shall exclusively deal with the (2+l)D-0(3) model 
where lattice sizes of the order of lO^x 10^ are computationally easy to handle and yet sufficiently 
large to observe essential features of the ordering process. 

Because the presence of topological textures implies the existence of additional scales rele- 
vant for the time evolution it is expected that phase ordering in such systems violates dynamical 
scaling. This has been demonstrated in previous simulations for 2D-0(3) non-linear and 
linear ^ sigma models. In two spatial dimensions the two- derivative (sigma model) term 
/ V^V^cPx is invariant with respect to the spatial scale, therefore in these minimal models 



the size of the defects is not fixed by the static energy functional. In fact, it was observed 
in |^,|^ that the corresponding length scales change with time and interfere with the correlation 
length scale that characterizes the aligning process. 

It appears desirable to include in the energy functional additional terms which stabilize 
the defects at a fixed finite size. Due to the scale invariance of the two-derivative term this 
requires at least two more terms which balance each other in the stable static configuration. 
The effect of a four-derivative term (which tends to increase the size of static structures) in the 
nonlinear (hard-spin) version of the model can be compensated by a (zero-derivative) Zeeman 
term which by itself tries to shrink local inhomogenities. Physically motivated by external 
magnetic fields coupled to the order field vector this term, however, explicitly breaks the 0(3)- 
symmetry and therefore prevents spontaneous alignment in random directions. It is only in 
the easy plane where the formation of disoriented domains can be observed as long as the 
field still has components in that plane. Simulations of phase ordering in such models have 
therefore mainly been concerned with the dependence of the resulting defect densities on the 
defect size 0. 

In the linear (soft spin) version of the model, however, where the length $ of the order- 
parameter field is not constrained, the (zero-derivative) potential V"($) can serve instead to 
set the length scale without breaking the 0(3) symmetry. With the familiar (<l>^ — /q)^ ansatz 
for V we arrive at a most simple model which allows for all features of spontaneous symme- 
try breaking, combined with the possibility of bag formation and the existence of localized 
structures with definite size embedded in the aligning field. Inclusion of a four- derivative term 
is necessary to prevent the collapse of these localized structures to zero size with subsequent 
unwinding. 

It is the aim of this work to investigate phase ordering in connection with the simultaneous 
formation of isolated topologically nontrivial structures stabilized by the energy functional. 
They appear as transient structures in every ordering process, but here they finally persist as 
stable 'particles' and 'antiparticles' with well-defined structure. We shall additionally consider 
the option that their net number is chosen as conserved observable. In section 2 the effective 
lagrangian is specified which comprises the minimal number of terms necessary to establish 
these features if we exclude all explicitly symmetry-breaking terms. In section 3 we briefiy 
discuss the stable static solutions of the corresponding energy functional which (in that model) 
only exist if the relevant coupling constant is below a critical value. Finally, in section 4 
we perform numerical simulations of individual 'events' which follow the ordering process in 
real time through a Langevin-type overdamped dynamics. For the exploratory purpose of 
this work we only consider the sudden quench scenario where a random initial configuration 
with correlation length less than the lattice constant is exposed to a low-temperature effective 
potential and a correspondingly small low-temperature stochastic force. 



2. The 2D-0{3) model with fourth-order stabiUzation 

We consider the 0(3)-symmetric lagrangian density in 2 + 1 dimensions in terms of the 
dimensionless 3-component field $ = $$ with $■$ = !, 

c=F^ Qa,$a-$ - T (^' - f^y - '' p^p") ■ (^) 

Apart from the usual sigma-model term this lagrangian contains the standard potential V{^) 
for the modulus field $ to monitor the spontaneous symmetry breaking, and a four-derivative 
('Skyrme') current-current coupling Pfj,p^ for the conserved topological current 
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pt' = —e'^'^p^ . (5^$ X dp^), (2) 



which satisfies dfj,p^ = 0. 

If we write the independent strengths A4 of the $^-coupling and C4 of the Skyrme coupling 
in terms of one common dimensionless parameter A and a length i 

A4 = X/f, C4 = Xf (3) 

then £ may be absorbed into the space-time coordinates. So, for A fixed, i sets the size of 
localized static solutions, and for continuous coordinates their total energy is independent of i. 
The overall energy scale is set by the parameter F"^. Of course, we are free to insert additional 
powers of the modulus field $ into the Skyrme term, the above choice being motivated to 
minimize interference with the $^ spontaneous symmetry-breaking mechanism. 

Having fixed the ^-dependence of the lagrangian as given in (|1]) and (0) we conveniently 
redefine the field and the parameters by 

^ = <^U\ F' = F^fl e = ef,-\ (4) 

This shows that for fixed i as /o goes to zero ( e.g. with increasing temperature ) the typical 
size i of static defects grows like I//0. This may be physically not unreasonable (cf. e.g. the 
discussion in the 3-dimensional case in |jlO[)- We omit the tildes in the following and absorb 
the €s into the length scale of space-time. Then we finally have for the static energy 

E = F'J (^d^^d,^ + ^ ($2 - 1)' + Xpopo) <fx. (5) 

For the following we will put the energy scale F'^ to unity. Note that the lagrangian (|l|) 
contains no symmetry-breaking term and in this sense is the close analogue to the massless 
chiral 3D-0(4) model. 

In the lattice implementation, we impose periodic boundary conditions for the field vectors 
which implies compactification of coordinate space to a torus S^ ^ S^. A stronger condition 
would be to require that $ is the same for all points on the lattice boundary, which would 
imply compactification of coordinate space to the two-sphere 5*^. In both cases the winding 
density po satisfies / PqcPx = B with integer winding number B. 
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FIG. 1. Dependence of the total energy on the coupling strength A for different values of B. 



3. Static soliton solutions 

Let us at first give a simple argument how for A less than a critical value the formation of 
localized bags will lead to field configurations with winding number B which are energetically 
more favorable than the standard Belavin-Polyakov (BP) soliton solution |TT| of the 0(3) 
nonlinear a-model where $ is confined to the 2-sphere $^ = 1 everywhere. For that purpose 
we consider idealized square-well bags, i.e. configurations with fixed total winding number B 
where $ is close to zero inside an area A and equals unity elsewhere, with all nonvanishing 
angular gradients confined to the inside of that area. For such configurations the winding 
density is p = Bj A and the first term in (^ does not contribute, therefore the bag energy (|^) 
is minimal for A = 2B and is obtained as 



XB. 



(6) 



Ebag 

For A > 47r this exceeds the energy E^ = AtiB of the standard BP solution , where the 
second term in (^) does not contribute. In that case the contribution of the Skyrme term can 
be scaled away by unlimited increase of the spatial scale. Therefore, for A > 47r idealized bags 
will not be stable but melt away into infinitly large BP-solitons. On the other hand, for A < 47r 
we may expect well defined stable bag structures with their spatial extent fixed by the choice 
of £, for any chosen value of B. 

The above argument for idealized square- well bags, involving only bulk energies, is indepen- 
dent of the shape of the idealized bag. For real bags due to the surface energy given by the first 
term in (^), the degeneracy of Ebag with respect to the shape will be lifted. With increasing 
values of B these surface effects will be less and less important. This is shown in fig.|l|, where 
the energy E{\) is plotted for 5 = 1 and B = 10. As expected, the corresponding energy 
curves E{X) lie slightly higher than the hnear result (H) derived above for idealized square-well 
bags, and approach that result with increasing values of B. 
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In numerical simulations on a discrete N x N lattice the lattice constant a defines an 
additional scale so we can expect independence of the energy from the scale i only as long as 
A^ ^ i/a = {ci/XiY^'^ ^ 1. Then the energy -E'(A4, C4) resulting from the static part of (|l|) will 
scale as E{\) in (|) of the single argument A = VA4C4, only. 

As £/a approaches 1 from above scaling violations set in. This is illustrated in fig.^ where 
Eb=i{X = 1) obtained on a square lattice with a = 1 is plotted for different values of i. Fig.^ 
shows that scaling holds with good accuracy for i > 4. 



FIG. 2. The ^-dependence of the energy Eb=i for fixed coupling constant A = 1. 

In the angular representation of the field 3-vector $ = $$ the length $ is the 'bag' field, 
and it is convenient to parametrize the angular part $ in terms of the profile function Q{x,y) 
and the azimuthal angle (j){x, y) with respect to some arbitrarily chosen cartesian basis 



($i,$2,$3) = $ (cos </) sin 9, sin (^ sin 6, COS O) . 



(7) 



For i = 10, i.e. well within the scaling region, the bag $(x,j/), the profile function Q{x,y), 
and the winding density p{x,y) of the resulting 5=1— configuration are plotted in figs.^]. One 
may recognize how the winding density is concentrated within the well-developed bag. The 
'profile'-function G drops from the value of vr in the center to zero outside the bag and the 
angular field (f){x,y) coincides with the BP- hedgehog form = arctan(?//x). So, although the 
angular configuration resembles closely the BP-soliton the energy Eb=i is (for A=l) only 3.70 
as compared to £^s=i = 47r. 
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FIG. 3. Profile function Q{x,y), bag ^{x,y), and the winding density p{x,y) of the 
B = 1-configuration. (A = 1) 
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FIG. 4. Total energy E = E^^ + Ehag, surface energy E^ = ^ J d-i^di^ cfx, bag energy 



Ebag = ^ J [\ (*^^ — 1) + PoPo) d'^x, and the minimum of the 'bag'-field <I>mini 
couphng constant A, for B = 4:. 



as functions of the 



The effects of the finite lattice constant on the field configurations can be studied as i 
approaches 1 from above: the bagfield $ develops a sharp dip by taking on a value very close to 
$ = only at one single lattice point while being a smooth extended function otherwise. The 
corresponding density p for a configuration with B = 1 then assumes the values of p = 0.25 
within each of the four adjacent lattice cells with p = everywhere else. Then, while the 
bagfield still is able to scale with i as ^e{x) = $(a;/£), the density p^ can no longer scale as 
pt{x) = p{x/i)/i'^. This then causes the scaling violations in the energy Eb=i as shown in 
fig.|^. If i is chosen still smaller, like i = 0.1, the density takes the value of 0.5 on two adjacent 
lattice cells while the bag zero disappears somewhere between the lattice points, such that the 
field $ «i 1 on all lattice vertices. In that case the lattice simulation produces a configuration 
which even looks as if the constraint $^ = 1 had been imposed. In that case large bags carrying 
multiple charges break up into individual B = 1 structures. In order to avoid such effects of 
the finite lattice constant we consider in the following scales i which are sufficiently large to 
be safely in the scaling region but still small enough for the resulting configurations to be well 
contained in a reasonably sized lattice (like N = 100 — 150). 

We finally proceed to the dependence of stable bag configurations and their energies on the 
coupling strength A. For B = 4 different contributions to the total energy and the depth of the 
bag profile are shown in fig.^], and in figs.|^ densities and bag profiles for i? = 10 and A = 1, 4, 10 
are compared. Referring to fig.^, one may still recognize nonvanishing bags for A > 47r as a 
consequence of the limitation of the lattice size: The size of the BP-soliton is restricted by the 
borders of the lattice, so that the contribution of the Skyrme term cannot really be scaled away 



for A > 47r. According to Hobart-Derrick's theorem the second and the third term in (|^) must 
contribute the same amount to the total energy, resulting in very flat but still nonvanishing 
bags for A > An. Figs.^ show that for B = 10 the bags already resemble the idealized bags 
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FIG. 5. Density and bag profiles for B = 10 and three different values of A = 1, 4, 10. {i = 10) 

discussed above quite closely, with a fiat interior in which the similarly fiat density is localized. 
Their radius is almost independent of the coupling constant A and fixed by the choice of i. 
The surface thickness depends on A but not in a dramatic way. Mainly the outermost tails of 
the bag profile are sensitive to A. However, increasing surface thickness squeezes the winding 
density towards the center of the bag such that the central value of the local density can be 
quite sensitive to A, especially for small values of A. 
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4. Relaxation after sudden quench 

Overdamped relaxation of initially random configurations leads to the formation of domains 
in which the field vectors are aligned in spontaneously chosen random directions. Size and 
orientation of these (dis)oriented domains change with progressing time, some of them growing 
on cost of others, such that altogether long-range order is increasing. Boundaries and edges of 
such domains are characterized by large angular field gradients. Energetically, large angular 
gradients favor formation of bags where the length of the field vectors deviates strongly from 
the vacuum value /q. Therefore the ordering process is accompanied by spontaneous formation 
of bags with winding density accumulated inside the bags. With progressing time these bags 
assume the spatial extent and profile dictated by the energy functional, while the areas of 
aligned field in which they are embedded finally grow and coalesce into a uniformly oriented 
vacuum. 

In order to follow this ordering process as it proceeds in real time we consider the equations 
of motion as obtained from (p, suppressing, however, all second-order time derivatives in 
comparison to a first-order time-derivative dissipative term: 

U = A$ - A ($2 _ i) ^ _ ^^2 d^pDid^ + ^^. (8) 

As long as we disregard second-order time-derivatives the damping constant 1/r which multi- 
plies the $ term can be chosen as unity, i.e. we identify the time unit with the relaxation time 
r. 

In (!]) we also have added a fiuctuational field $,{x, t) to represent gaussian white noise; in 
principle its presence and strength is dictated by the dissipation-fluctuation theorem. Inde- 
pendent from this stochastic dissipation other damping mechanisms could be present, like the 
rapid cooling due to the Bjorken expansion of a hot hadronic fireball |]T^. So, depending on 
the specific physical situation the damping rate and other parameters in (P) may be subject 
to an appropriate time (or temperature) dependence. Here, however, for definiteness, we will 
keep them fixed during each individual evolution. This corresponds to a sudden quench where 
the initially (at t = 0) hot configuration is exposed for t > to the low-temperature (T = 0) 
effective action. Consequently, we also generally will omit the noise term ^. Average results 
are not sensitive to it, anyway; only the accidental features of late-time configurations reached 
in individual evolution events are affected by the noise term. 

Initial configurations are chosen such that at each lattice vertex (i,j) {i.,j = 0,..,N) the 
field vectors $ point in some random direction i.e. at each point of the lattice the angle (f){i,j) 
is selected randomly from the intervall [0, 2tt], the angle 0{i,j) from the intervall [0, vr]. Through 
this choice the finite lattice constant acquires physical meaning as providing a measure for the 
magnitude of the initial correlation length. The moduli $(i, j) of the field vectors $ are chosen 
as absolute values of a Gaussian deviate around the symmetry center $ = 0. Again, late-time 
average features of the resulting configurations do not depend significantly on the mean square 
deviation of this initial gaussian distribution. In fact, similar results are obtained even if the 
initial configuration is constrained to the 2-sphere $^ = /q. The reason for this is that the 
system reacts to the initially large local angular gradients by reducing the length of the field 
vectors almost everywhere to values which are small as compared to /q. This happens early 
during the first few time steps, accompanied by some next-neighbor alignment. Therefore the 
initial length distribution is almost instantly forgotten. 
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At the borders of the lattice periodic boundary conditions are enforced. If we divide each 
elementary lattice cell with lower left corner {i,j) into two triangles (e.g. by the same diagonal 
in all cells), then the map ^{i,j) maps each triangle onto a spherical triangle on the sphere 
$^ = 1 cut out of the surface of this sphere by the (shortest) geodesies which connect the image 
points of the corners of each triangle. The local winding density p{i,j) then is defined as the 
sum of the (oriented) areas (divided by Att) of the two spherical triangles which form the images 
of the lattice cell with lower left corner {i,j). Because the area of each spherical triangle is less 
than 27r each square lattice cell can contain at most one unit of total winding number. The 
periodic boundary conditions guarantee that the total winding number 

N-l 

B=Y.Pi^^j) (9) 

i,j=0 

summed over the whole lattice is integer, so that initial configurations can be selected with 
some desired integer value of B. We also define the "number of defects" 

N-l 
i,j=0 

by summing up the absolute values of the local winding densities. Of course, for random or 
slowly varying smooth configurations D generally is not an integer, but if a configuration de- 
scribes a distribution of localized defects (and antidefects) which are sufficiently well separated 
from each other, then D is close to an integer and counts the number of these defects (plus 
antidefects). In that case we can define the numbers A^+, A^_ of "particles" and " antiparticles" 
through 

B = N+-N_ D = N+ + N_ . (11) 

We shall, however, in the following (sloppily) call D the "particle number" (even if it is not 
integer). (Alternatively, D could be defined as the sum of the absolute values of the areas of 
all spherical triangles considered above. For a random configurations this would result in an 
average value of {D) = iV^/4 (Kibble limit [|T^). Our definition (|TD|) for random configurations 
leads to (D) = 0.73A^^/4 which implies a slightly different definition of the initial correlation 
length.) 

The total winding number B always is integer and occasionally will undergo discrete jumps 
in the update sweeps. For well-developed localized structures this corresponds to unwinding 
defects or antidefects independently, such that B decreases or increases by one or more units. 
This 5- violating propagation is characteristic for the trivial topology of the linear 0(3)-model. 
However, with the evaluation of B for each instantaneous configuration we may in the lattice 
simulation implement an (optional) B filter which in each time step rejects configurations that 
violate B conservation. This eliminates all independent unwinding processes. Only simulta- 
neous annihilation of defect and antidefect in the same time step remains possible, and, as 
the update proceeds locally at each lattice vertex it can happen only if defect and antidefect 
overlap. This 5-conserving evolution is characteristic for the nontrivial topology of the nonlin- 
ear 0(3)-model. Of course, we expect that severe differences between both types of evolutions 
appear only if D is comparable to B. 
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FIG. 6. Time evolution of the correlation function for B = and A = 5. The time is given in units 
of the relaxation time r. 



In order to produce configurations with well-developed stable bags we choose for the fol- 
lowing quench simulations values for the coupling constant between 0.5 < A < 5. This is well 
below the critical value of A = 47r. 

For most of the results presented below the scale i is chosen as £ = VTO; this is small enough 
to allow for the formation of numerous bags on a reasonably sized lattice (like 150 x 150), but 
is still close to the onset of the scaling region (cf. fig.0) to suppress effects of the finite lattice 
constant. 

To obtain a more quantitative measure for the size of ordered domains we consider the 
correlation function 



N N 
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ij=0 k,l=0 
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i j=0 k,l=0 
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where the k, I sum is restricted such that the distance r = J{k — iy + (/ — j)^ between lattice 
vertices (z, j) and {k, I) lies inside bins of unit size around fixed positive integers R. The typical 
shapes of these correlation functions are shown in fig.^ for an evolution with ^-conservation 
on a 150x150 mesh, for increasing time. For the initial configuration C{R) vanishes for all 
R > 1 which reflects our choice of the initial correlation length. In the very early part of the 
relaxation process [t < 1) these correlation functions approach zero within less than 5 lattice 
units and stay close to zero for larger distances. For t > 10 they drop to small values above 
10 lattice units, but stay positive with small oscillations ( with wavelengths of more than 20 
lattice units, increasing with t). For t > 200 these oscillations no longer fit into the 150x150 
lattice, the correlation functions show a monotonous decrease (for R < N/2). Their minima 
(near R = N/2) increase towards unity, which indicates that finally the extension of the ordered 
domains reaches the size of the lattice. 
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FIG. 7. Total energy E, 'particle number' D, 'radius' Rd of ordered domains, and the length ($) 
of the field vector averaged over the whole lattice, for a few B = 0-conserving evolutions (for A = 1, 
i = vlO, on a 120 x 120 lattice). The time is in units of the relaxation time r. 



Comparing with typical field configurations during the relaxation, the half- maximum dis- 
tance, i.e. the distance R where C{R) drops below 0.5 apparently provides an appropriate 
measure for the 'radius' Rd of ordered domains. Of course, this is a rather arbitrary and not 
very precise convention, but it captures the essentials of the ordering process in view of the fact 
that the boundaries of the ordered domains are not sharply defined. A typical feature of these 
correlation functions is the appearance of a shoulder for small distances {R < 5) for late times. 
This reflects the formation of the ordered textures, i.e. spatially extended angular twists which 
locally prevent alignment of the field vectors over distances of the order of i. Naturally, this 
effect gets especially prominent for evolutions which proceed in configurations constrained to a 
large total winding number B (see below). 

For a series of 5-conserving evolutions which start off from different randomly chosen initial 
configurations, selected however for winding number B = 0, fig.0 shows the time dependence 
of the 'radius' Ro of ordered domains, the total energy E, the 'particle number' D and the 
length ($) of the field vector averaged over the whole lattice. 

Figs.^ show the typical features of the field in an 70 x 70 section of a 150 x 150 lattice during 
different stages of such an evolution after t = 1, t = 10, and t = 1000 relaxation time units. 
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FIG. 8. Typical features of the field during different stages of an evolution (for A = 5, with B 
conserved at B=0). In the left column the lengths $ are plotted over the spatial x-y plane which 
provides a 3d-view of the momentaneous bag structures. The right column presents the projection of 
the ^-vectors at each lattice site on the <l>i — <1>2 plane. 
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FIG. 9. Growth of the size Ro for different values of the couphng constant A (A 
B = 0-conserving evolutions on a 120 x 120 lattice. 



0.5,1,2,5) for 



Fig.|^ shows the growth of the size Rd for different values of the coupling constant A (A = 
0.5, 1, 2, 5) for B = 0-conserving evolutions which start off from identical initial configurations. 

One can distinguish three phases of the ordering process: 

i) During an initial 'relaxation' period which takes a few (relaxation-) time units, the lengths 
$ of the field vectors initially rapidly decrease and then vary around small values of about O.l/o, 
the number of defects D drops from its starting value (which is of the order of 0.73A^^/4) by 
about one order of magnitude, accompanied by a corresponding loss in total energy. During 
this period the growth of Rd closely follows a power law 



R 



D 



at", with a ^0.4 ±0.01. 



(13) 



This exponent is with good accuracy independent of the coupling constant A, of the initial 
configuration, and of the lattice size N (as long as A^ ^ 1). It is also independent of the 
scale parameter L By the end of this period ordered domains extend over several lattice units 
{Rd ~ 5). 

ii) The second phase (10 < t < 100, depending on A and £) could be termed the 'roll-down' 
phase. It is characterized by the increase of the (spatial-)average length ($) of the field vectors 
towards the vacuum value /q. Actually, this roll-down process is rather slow; it takes of the 
order of ~ 100 relaxation time units for the space-averaged ($) to approach /q. Locally, this 
increase of $ happens only in the interior of ordered domains, which results in the formation of 
numerous dense and initially often connected bag structures located around the boundaries of 
these domains. So, during this phase it is evidently the $^ potential which drives the evolution. 
Therefore the onset of this second phase depends on the scale parameter £ and on the coupling 
strength A. Large values of £ prevent the onset of this bag forming process for a long time, so 
the angular alignment proceeds further while ($) is still small. This results in much smaller 
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total particle numbers when the bags finally are formed. Similarly, this second phase starts 
earlier for larger values of A. 

Interestingly, the growth rate of the size of the ordered domains remains basically unaffected 
by this roll-down of $: the increase of the correlation functions proceeds monotonously through 
this phase. There is, however, an effect on R£, from the shoulder which appears in C{R) for 
5 < i? < 10 due to the developing localized extended winding structures. If this shoulder passes 
through the halfmaximum which is used to define Rd it leads to a deviation from the power 
law (|T^) which is especially pronounced if the bag formation sets in late, (i.e. for small values 
of A or large £), when the size of the domains in which $ approaches /o is larger. When the 
bags are fully developed the ordered-domain size has increased by about a factor of two, so the 
bags then are embedded in a patchwork of ordered domains (with $ k. /q) which extend over 
10-20 lattice units. 

iii) The further development proceeds by the bags slowly moving around, eating up smaller 
ones or uniting with others they meet on their way, or, annihilating with others of opposite 
winding number. They assume sizes which correspond to the chosen scale £ and refiect the 
partial winding number contained in their interior. So, naturally, the evolution during this 
late period depends more on accidental features of the individual configurations as they have 
developed up to that point. However, it is interesting to note, that on the average, the growth 
of the ordered domains (now with $ = /o) again approximately follows the power law ( [TB| ) 
with a around 0.4 (see fig.|^). Finally, if Rd has reached values near or geater than N/A the 
finite size of the lattice (with its periodic boundary conditions) affects the long-range part of 
the correlation functions resulting in a rapid artificial increase of Ro. 

The optional filter on the total winding number B allows to compare evolutions with B 
conserved at some initial value with others where B may jump freely during the course of 
the relaxation. As long as the actual values which B takes on are small as compared to D 
there is almost no difference between 5-violating and 5-conserving evolutions. For the major 
part of the evolution the local particle-plus-antiparticle density is so high that the evolution 
is dominated by annihilation processes, and occasional unwinding jumps in B by one or two 
units play no significant role. Only at late times (t > 500) when D has dropped below a few 
percent of N"^ the possibility of spontaneous unwinding makes a noticeable difference and it is 
accompanied by a correspondingly more rapid increase of the size Rd of ordered domains. 

On the other hand, by choosing a large initial value for 5, the S-conserving relaxation 
allows to study the formation of ordered domains in a 'baryon-rich' environment. We present 



a series of such events in fig.|TO| with B > 100 on an A^ = 120 lattice (for A = 1). In the 
early stages as long as D ^ B there is almost no difference as compared to the B = case. 
However, for t > 10 D approaches the value of B fixed at B > 100. At these times the value 
of $ is still much smaller than /q at most lattice points and the constraint on winding number 
presents a severe obstacle for the growth of many aligned domains, i.e. it causes a strong 
deviation from the power-law (|13]). In fact, the ordering proceeds in such a way that only very 
few oriented domains start growing with $ approaching /o in their interior. These few aligning 
domains squeeze the regions with nonvanishing winding (and very small $) into coherent large 
bags, such that the whole space becomes separated into large aligned areas and large bags. 
This structure is refiected in the correlation functions as a very pronounced shoulder in the 
range around 5 < i? < 10 which rises with increasing time due to the further increase of long 



range correlations (cf. fig. 10). Evidently, near t ~ 100 this leads to an almost instantaneous 
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FIG. 10. Time evolution of the correlation function C{R), and a series of evolutions as in figj^ 
however in a 'baryon-rich' environment where B is selected and conserved at values B > 100. 



strong increase in Ra, which thus appears more as a consequence of the definition of Rd as the 
half-maximum distance rather than an actual abrupt increase of the size of aligned areas. The 
further development then is characterized by the formation of one large bag which comprises 
almost all of the winding number (cf. fig.pT|) , while alignment in the surrounding Vacuum' in 
most cases progresses slowly according to (|13D with a ~ 0.4 before finite-lattice-size effects set 
in for t ~ 500. 
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FIG. 11. Typical example of the field in a baryon-rich environment (for A 
after t = 1000 relaxation time units. 



1, with B > 100) 



Apart from fixing the spatial extent of the finally formed bags the scale parameter £ affects 
during the early stages of the evolution the duration of the first period where the bags are not 
yet fully developed. Thus it allows to monitor the total particle number present at the time 
of bag formation. This is an interesting aspect for evolutions where one does not consider a 
sudden quench but allows /o to change with the temperature of the system. According to eq. @ 
this transforms into changing i with time. 
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5. Conclusion 

We have presented here numerical simulations of phase ordering for a 3- vector field through 
spontaneous symmetry breaking in two spatial dimensions, with specific attention to the inter- 
play between the aligning process and the formation of stable extended topological structures. 
For that purpose the effective action is chosen in such a way that the localized topological 
defects which necessarily accompany the formation of randomly oriented aligned domains are 
stabilized with a definite size as ordered localized structures embedded in the aligning field. 
The stabilizing terms do not break the 0(3) symmetry explicitly. For definiteness we have 
only considered the sudden quench scenario and assumed overdamped dynamics. Apart from 
the interesting features of the coarsening transition such processes may serve as models for the 
spontaneous creation of extended particles and antiparticles or clusters of those out of a hot 
random field ensemble. 

Three aspects are of peculiar interest: The growth rate of the size of aligned domains follows 
a power law with an exponent of approximately 0.4, which persists through the early relaxation 
phase and the subsequent roll-down phase and is apparently quite independent of the stabilizing 
terms in the action. This exponent is in agreement with previous results found in models 
without (or with only one of the) stabilizing terms for the average defect-defect separation [Q 
and the spin-spin correlation . The formation of stable defects is most prominently reflected 
in the shape of the equal-time angular correlation functions. During these earlier phases the 
ordering process is dominated by defect-antidefect annihilation which reduce the initial particle 
number by up to two orders of magnitude before the remaining defects slowly take on their 
stable conformation. By that time (which is of the order of several hundred relaxation times) 
the final alignment of the remaining few large disoriented domains depends sensitively on the 
accidental spatial configuration of the few surviving extended particle clusters. 

The scale parameter £ which determines the size of the resulting defects does not affect the 
early relaxation period but it has a pronounced influence on the subsequent evolution: larger 
values of £ suppress the increase of $ towards its vacuum value /o over larger spatial areas and 
thus delay the onset of the roll-down phase. This leads to a reduction of the final density of 
particles plus antiparticles when the bags finally emerge. Through (^ C. is directly related to 
/o which implies sensitivity of the final total particle number to the quench velocity. 

The third interesting feature concerns 5-conserving evolutions in an environment with large 
values of B. Similar to the case of large i the onset of the roll-down is delayed until most of 
the possible annihilations have taken place, i.e. until the total particle number approaches the 
fixed winding number. The subsequent growth of only a few aligned domains squeezes regions 
with nonvanishing winding density into large coherent particle clusters which fill the interior of 
large bags. 

Extending the present considerations to the 3D-0(4) model appears as a challenging task 
in view of the ongoing discussion of the chiral phase transition, the formation of disoriented 
chiral domains, and baryon-antibaryon production in the cooling of hot hadronic plasma. It 
may help to establish further links between assumed quench scenarios and signatures to be 
expected from emitted particles. 
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